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Abstract 

A system of a metastable phase with several sorts of heterogeneous centers is consid- 
ered. An analytical theory for the process of decay in such a system has been constructed. 
The free energy of formation of a critical embryo is assumed to be known in the macro- 
scopic approach. At first all asymptotic cases are investigated and then a general interme- 
diate solution is suggested. All approximate transformations are accomplished with the 
corresponding numerical estimates and analytical justifications. This is the second part 
of the manuscript, the first part can be also found in this archive. 
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This is the second part of the manuscript, the first part can be found in the cond-mat e-print 
archive at the same publishing date. 

1 Intermediate approximate solution for several sorts 
of centers 

In the case of nucleation on one sort of centers the expression for g%(z) is unknown only when 
the spectrum is cut off by exhaustion of heterogeneous centers. So, one doesn't know it in 
the situation when the converging force of the operator Si is very strong. In the case of the 
several sorts of centers the situation is another. One doesn't know every term in Y^j9j( z )- 
Then one can come to the situation where according to the first iteration the spectrum is cut 
off by the exhaustion of supersaturation initiated by droplets formed on some sort of centers 
but in reality all heterogeneous centers of that sort are exhausted and those droplets practically 
don't consume the vapor (due to their quantity). So, it is necessary to have some more precise 
expression for gi which allows to calculate the next iteration for #j. 

1.1 Mono-disperse approximation 

The length corresponding to the cut-off of the spectrum by the substance exhaustion is practi- 
cally one and the same for all sorts of the droplets (all Tj have approximately one and the same 
valued) . Analyzing the subintergal expression in the equation for in the pseudo homogeneous 
case one can see that the subintegral expression s = (z — x) 3 exp(— I 1 * J2j 9j/ $*) connected with 
a variation of the supersaturation in the formula for is a very sharp function of x. It is less 
than the function 

Sbel = 9(z - X){Z -Xf (1) 

and greater than the function 

s ab = Q{z-x){z-xfeM- T -^^ X ^) ■ (2) 

Fig. 6 illustrates the behavior of the values s a b and Sbei- It can be seen that they practically 
coincide. This is the reason of applicability of the first iteration as a good approximation for 
the precise solution for one sort of centers. 

Let us introduce an approximation for this function. One has to extract the region of the 
sizes of the droplets which are essential in the vapor consumption. This consumption is essential 
when 

z « Ax (3) 

where Ax is the length of the cut-off by the supersaturation. Certainly, the sizes of the droplets 
which are essential in the vapor consumption must be smaller than Ax because the effectiveness 
of the iteration procedure in the homogeneous case is based on the fact that the droplets formed 
at the almost ideal supersaturation determine the formation of the spectrum. One can see that 
s is a very sharp function of x and this is the reason for the monodisperce approximaion for s. 
Now one has to define characteristics of this approximation. 

1 Alternative situations correspond to the asymptotic cases and have been already described. 
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Figure 6. 
Functions s ao and s^i. 



The condition for the differential width is that the subintegral function falls two times in 
comparison with an amplitude value. For the differential half-width 81/2X one can get the 
following expression 

S 1/2 x = (f - -rj^)x . (4) 
The integral half- width Ay 2 x can be obtained from the corresponding equation 

N ess x J* %—^n^ , (5) 

where iV ess is the characteristic number of the droplets obtained as iV ess = /* iAx^xn^ and 
the r.h.s. of the previous equation is the first iteration approximation. This gives 

A 1/2 x = -x (6) 

and it practically coincides with 5i/ 2 x. 

The subintegral function s is now decomposed into the essential part where 

z 

and the tail where 

z 
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The tail will be neglected and due to the small relative size of the essential part the mono- 
disperse approximation of sizes for the droplets in the essential part will be used. As the result 
one can get an approximation 

g(z) = ^-^z 3 , (7) 



where N(z/4) is the number of the droplets appeared from x = till x — z/4. 

Since the spectrum is cut off by the exhaustion of the supersaturation in a frontal (sharp) 
manner, the value of gi is unessential before z = Ax because it is small. After the moment of 
the cut-off it is unessential also because there is no formation of the droplets. So, instead of 
the previous approximation one can use 

N i (A l x/4) 3 

9i(z) = z A . (8) 

Now let us turn to the heterogeneous case. The exhaustion of heterogeneous centers makes 
the subintegral function more sharp and the mono-disperse approximation becomes at x ~ AjX 
even better than in the pseudo homogeneous situation. Certainly, one has to use N(Aix/4) 
calculated with account of exhaustion of heterogeneous centers (but at the coordinate, obtained 
without any account of exhaustion of heterogeneous centers). 

Fig. 7 illustrates the form of the spectrum on the base of "mono-disperse approximation". 
The case of the pseudo homogeneous condensation is considered. Three curves are drawn: the 
spectrum in the mono-disperse approximation f appr ; the spectrum in the first iteration f\ which 
can be considered as a very precise approximation and the numerical solution /. Two of them 
/ and fi coincide and correspond to a thick line. The back side in the first iteration is more 
sharp than in the "mono-disperse approximation". Nevertheless the deviation isn't so essential. 
It can be eliminated by a perturbation theory. 

One can get N(Aix/4) from the solution of the equations for the separate condensation 
process because only the lowest length of the cut-off is nesessary. This length is given without 
any cross influence taking into account due to the frontal character of the back side of the 
spectrum. Then one has to repeat this procedure with the already known iVj. 

1.2 Final iterations 

Now one can see a way to solve the system of the condensation equations. 
A. Method of separate lengths. 

At first one has to solve equations for all separate nucleation processes 

9% = f*i[ (z- xf ex.p(-r~)6idx = Gi(g i} 0;) , (9) 
Jo <P* 

0i = exp(-/* / exp(-T~)dx) = Si(gi) (10) 

Vtot i Jo 

for every %. This solution is given by the standard iteration procedure. In the second approxi- 
mation one can calculate the value of 6{2){Aix/4) explicitly 

e im {A lX /4) = eM-f* i— / AW4 exp(-£j^ 4 )^) = exp(-/, ^{^)- 1/A C) , 

Vtot iJ0 4<P* Tjtot i 4$* 
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Functions f, fi and /< 



where 

rl/4 

C= / exp(-z 4 )dx « 0.25 . (12) 
Jo 

Then 

iV l(2) (A^/4)=^ ,(1-^(^/4)) . (13) 

One has to fulfil these calculations for every sort of the heterogeneous centers. Considering 
these approximations as the initial ones it is necessary to do only one more step of the iteration 
procedure to get the suitable results. They will be marked as the 'final' ones. One has to 
calculate 0j final (oo) and N fina i i(oo) 

9> fmal (z) = exp[-/, ^ f^M- ^^^^^ ^dx} , (14) 

N fmal ,(oo) = THot ,[1 " exp(-/, ^( Ejr ^ (2 i (AW4 V 1/3 ^] • (15) 

These expressions are valid under the reasonable separation of heterogeneous centers into 
sorts when the centers with approximately one and the same height of the activation barrier 
are considered as one sort. Meanwhile it is obvious that when one splits one sort into very 
many sub-sorts then one can formally attain a wrong result. 

Let us suppose that the number of heterogeneous centers ensures hi ^> 1, i.e. the fall of the 
supersaturation leads to the interruption of formation of the droplets. One can split this sort 
into so many sub-sorts (j) that for every sub-sort hij <^ 1, i.e. one can see the exhaustion of 
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the sub-sort in the separate process of the condensation. Moreover, one can assume that this 
exhaustion is finished before z attains the quarter of the length of the spectrum initiated by 
the fall of the supersaturation. Since A^x is proportional to ^Yot (where i]^ tot is the total 
number of the centers of the given sub-sort) and SijX doesn't depend on this quantity, the 
required property is obvious. Then after the summation over all sub-sorts one can see that 
the total number of the droplets formed up to Aix/4 coincides with the total number of the 
heterogeneous centers. This conclusion is wrong. 

The reason of the error is that the width of the spectrum is smaller than the width in the 
process of the separate formation. Evidently, the mono-disperse approximation doesn't work 
at such distances. Under the reasonable definition of the heterogeneous sorts all characteristic 
lengths are different. So, ordinary, there is no such effect in this situation. But still one has 
to correct these results. The simplest way is to consider all centers with the difference in the 
activation barrier height less than one thermal unit as one and the same sort. This solves the 
problem, but this requirement is rather artificial. 

Now these results will be corrected in more elegant way. The mono-disperse approximation 
leads to some already defined functional form for g and for The form is already known and 
now only parameters in these functional dependencies have to be determined. 

B. Method of the unique length. 

The unique length of the spectrum Ax will be chosen. At this very moment Ax is unknown, 
but satisfies the following inequality: 

Ax < AiX 

for every sort of the heterogeneous centers. This leads to 

9i(z)= 1 ' V , (16) 

which gives 

N t{2) (Ax/A)=r ]ltot (l-e m (Ax/A)) , (17) 
Oi (2 )(Ax/4) = exp[-/, / exp(-- £ N^Ax/A^dx] . (18) 

Vtot i JO ( i'*^oo j 

The weak dependence Fj ~ T on the sort of the centers is neglected here for simplicity, while 
the strong dependence /*j on the sort of centers is taken into account. After the substitution 
one can get the system of the algebraic equations for N^Ax/A) 

rr, rAx/4 F r 3 

N, t (Ax/A) = th tot (l - exp(-/, / exp(-- £ N 3 {Ax / A))dx)) . (19) 



Let us simplify the last system. To calculate the integral one can note that 

rx 

/ exp(-x 3 )dx w x x < 1/4 . (20) 
Jo 

So, the trivial dependence of the r.h.s. on Nj disappears. One can come to 

iV,(Aa;/4) = Vi tot (l - exp(-/, *-^-^)) . (21) 

Vtot i 4 
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On the other hand one can use the sense of Ax as the half-width of the spectrum due to 
the fall of the supersaturation: 

At 3 r 

— 5>(A*/4)- = l . (22) 



After the substitution one can get the equation for Ax 

^£Wl-exp(-/*^^))^ = l • (23) 

"oo j '/tot j ^ 



To solve this algebraic equation one can choose the initial approximation as a solution of 
the following equation 

At 4 T 

E/.isr = i ■ ( 24 ) 



4 ^ J <£> 

j * 

Iterations are constructed according to 

w(l -exp(-/, °° (t 1 ))- = 1 . (25) 

'too j '/tot 3 ^ ^* 

When the last equation is solved then equation fl2~T|) expresses iVj (Arc/4) through Ax and 
solves the problem. Equation (16) gives the expression for g^ and for the supersaturation as 
the function of time. 

One can estimate the relative error by the pseudo homogeneous case which leads to 

I iVj(oo) - N f ina , Aoo) I I A - B I . . 

t™**\_LL<l L „ o.02 (26) 



AT*(oc) 

The investigation of the process of nucleation on several sorts of heterogeneous centers is com- 
pleted. 



2 Solution for heterogeneous centers with quasi con- 
tinuous activity 

2.1 Iteration procedure 

The smoothness of r](w) is supposed here. It allows to hope that the cross influence of exhaustion 
of one type of heterogeneous centers on the cut-off of the droplets formation on another sort of 
centers will be not so strong and one can apply some modifications of the standard iteration 
procedure. It can be constructed in a following way: for initial approximations one can choose 

g o {z,w) = , O = 1 ; (27) 
the recurrent procedure is defined according to 

g i+1 (z,w) = G w (gt ot ,0i(w)) , gf\z) = J dw 9i (z,w) . (28) 
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e l+1 (z,w) = sugf) (29) 

The monotonic properties of G w , J dw and S w lead to the chains of inequalities analogous 
to (??)-(??), which ensure estimates for precision of approximations. Certainly, to calculate 
iterations one has to use some expression for r) tot (w). So, at first one has to consider the limit 
cases and obtain some estimates. 



One can see the absence of r\ tot in the r.h.s. of (|29|) . So, for arbitrary r) to t the power 
of exhaustion will be determined only by w. The centers with the high activity w are almost 
exhausted during the process of the condensation. The centers with the relatively small activity 
w remain free. The intermediate region has rather small size. To give the qualitative estimates 
one can see that the supersaturation ( appears only in the function f^. For one can obtain 
the following estimate 

ft ~ /* exp(— constx 6 ) 3 < e < 4 (30) 

This estimate goes from the obvious fact that the spectrum of the droplets must be wider than 
the mono-disperse one and the intensity of nucleation must decrease in time. On the base of 
the last estimates one can see that the width Ax of the spectrum is a well defined value. 

The rough approximation for exhaustion of heterogeneous centers can be obtained from the 
first approximation in the iteration procedure 

9 = exp[— /* exp(\w)n 00 z] (31) 

and for the value at the end of nucleation period 

final = exp[-/„exp(Aw)n 00 Ax] . (32) 
One can define w < and w > according to 

w< = w - - , w> = w + - , (33) 

where 

1 . In 2 . 

w = T-H- f T-) , e~l ■ (34) 

These constructions are illustrated in Fig. 8. For w > w > almost all heterogeneous centers are 
exhausted. For w < w< almost all heterogeneous centers remain free. This remark allows to 
rewrite the expression for g tot in the following form 

/■u?> roc 

g tot = dwg(z,w)+ r] tot (w)dw . (35) 

The size of the intermediate region corresponds to a variation of w of an order 1/A or to the 
variation of AF in one thermal unit. Since AF ^> 1, the relative variation of AF in the 
intermediate region is small. So, it is reasonable to put in this region 

r)tot(w) = V* = const ■ (36) 

One can spread this approximation over the region w > w > because the behavior of r\ tot for 
these centers isn't important. Only J f]totdw over this region is essential. Certainly, one has to 
get the boundary w max of this region by 

rWmax POO 

^dw = r]*(w max - w>) = / r] tot (w)dw . (37) 



s 
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Figure 8. 

Function 9 and boundaries w < , u>>. 



For w < w < the accuracy of this approximation is not so essential because these centers remain 
free. 

After the scale transformations Tg tot /$* — > G; Tg /$* — > g 



9 = 



exp(Aw) P \z - xf exp{-G)9dx , (38) 
•P* Jo 

/ 1^771 ax 
dwg , (39) 
-oo 

9 = exp[— f*noo exp(Xw) / exp(— G)dx] . (40) 

Jo 

By the appropriate shift of w one can put w max to the zero. In the appropriate scale of w one 
can put A = 1. In the appropriate scale of z one puts the coefficient in ( pHf ) to unity. As the 
result, only one parameter - the coefficient /*^ 4 exp(3Xw max /4:)^V 4: / (T 1 / 4 ^ 4 ) remains in the 
last equation. One can mark it by A and get 

g = exp(w) / (z - x) 3 exp(-G)9dx , (41) 
Jo 

G = f dwg , (42) 

9 = cxp\-Aexp(w) [ exp(-G)dx] . (43) 
Jo 
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An iteration procedure can be constructed as 

g i+ i = exp(w) I (z- x) 3 exp(-G i )6'jdx , (44) 
jo 

Gi = / dwgi , (45) 



oo 



9 i+ \ = exp(— Aexp(w) / exp(— Gi)dx) , (46) 

Jo 

^o = , Q = 1 . (47) 



In the first iteration 



(?! = exp(w)— , (48 



G?x = f efc^i = i- , (49) 

J —oo 4 

6 1 ! = exp(— Aexp(w)z) . (50) 
The second approximation gives for 6 the following result 

02 = exp[— Aexp(w) / exp( )dz] , (51) 

jo 4 

then for the final value 

fin ai 2 = expt-AexpHl.28 1 / 4 ] (52) 

and for N tot 

N tot = [ [1 -exp(-Aexp(w)l.28 1/4 )]dwri* . (53) 
For g in the second approximation 

g> 2 = exp(w) / (z — x) 3 exp(— Aexp(w)x) exp(— x A /A)dx . (54) 
Jo 

The analytical expression for gi and further iterations can not be calculated. 



2.2 Universal solution 

The system of equations ( f4i"D - ( [43]) doesn't allow the universal solution as in the case of 
homogeneous condensation when ([4l"D hasn't 9 in the r.s.h.. But the operator in the r.h.s. 
of ( [4"1"D ensures a rapid convergence of the iterations. The worst situation for the iteration 
convergence is A = 0. In this situation one has after rescaling the universal system 

g = exp(w) / [z — xf exp(— G)dx , (55) 



G 



f° 

/ gdw . (56) 

J — oo 
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Figure 9. 

Universal form of spectrum in the pseudo homogeneous case. 



This system has the universal solutioni which will be marked by Go- This universal form of 
the size spectrum is shown in Fig. 9. 

The evolution is determined by the first three momentums (and the zero one) of the distri- 
bution function 

Hi(z)= f x* exp(-G)6dx , 2 = 0,1,2,3 (57) 



and in the pseudo homogeneous case it is determined by 

fii(z) — / x l exp(— G )dx . (5£ 



After the end of the short period of intensive formation of droplets one can substitute in 
fl57D and ( |5Tj| ) oo instead of z in the region of the integration. So, the further evolution will be 
determined by four constants! /i«(oo) i = 0, 1, 2, 3. If one takes for G the universal solution Go 
(when A = 0) then the values of /ij are the universal constants which can be obtained by the 
unique numerical solution of the last system. 

Now one has to return to the iteration procedure. One can use Go as the initial approxi- 
mation in the iteration procedure and get 

6>! = exp[— Aexp(w) / exp(— G )dx] , (59) 
Jo 

#2 = exp(u>) / (z — x) 3 exp(— Go) exp(— Aexp(w) / exp(— Go)dx')dx , (60) 
Jo Jo 



2 The trivial integration over w leads to the equation analogous to the homogeneous case. 
3 The main influence is due to /j,q. 
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/O rz rx 

exp(u>) / (z — x) 3 exp(— Go) exp(— Aexp(w) / exp(— Go)dx')dxdw . (61) 
-oo Jo Jo 

The decomposition of exponents gives 

oo (_A)i 

«-g^<W , (62) 

where z 

Pi= [ {z - xfexp(-G )jUx)dx , (63) 
Jo 



Joo(x) = / exp(-G )dx . (64) 
Jo 

Then the function B(z) = J Z exp(— G2{x))dx at z = oo will be the following 

POO POO 00 ( /4V 

B(z = oo) = A exp(-G 2 )dx = A / exp(-P (z)) exp(- V ^ '—PAz^dz . (65) 

JO Jo ^ (2 + 1)! 

The decomposition of the last exponent leads to 

00 00 ( — 1)3 (—A) ij 

where 

(7« = / exp(-P (x))P t J dx (67) 
■/ o 

are the universal constants. For the total number of the droplets with a given activity 

N to t{w) = ?7*(1 - exp(-5(2; = oo) exp(u>))) . (68) 

In the same manner one can take into account the deviation of 77* from a constant value. 
After the decomposition of 77* into Tailor series one can take the initial approximation as Go 
and act in the manner similar to (|55| ) - (p^). As a result, some expansions in the powers of the 
parameter A, the derivatives of r\tot and the universal constants can be seen. 

Fig. 10 illustrates these constructions. Here the most interesting function I 00 = J °° exp(— G)dx 
is drawn. Two curves are shown, the lower curve corresponds to the precise numerical solu- 
tion, i.e. to J °° exp(— G)dx. The upper curve corresponds to B{z = oo)/v4 = J Z exp(— G 2 )dx. 
This curve is very close to numerical solution. It is very important that both curves resemble 
a straight line. So, already the linear approximation on A corresponding to the account of 
the first leading (linear) term will be good. It means that in the last expression one can take 
instead of B{z = 00) the value A f °° exp(— Go)dx (on the base of the first iteration one puts 
~ A Io°° fidx). 



2.3 Wide spectrum of active centers 

Three regions were extracted in the spectrum of activities: the region of active centers (they 
are almost exhausted), the region of centers with small activity (they remain practically free 
during the whole process) and the intermediate region. In the majority of the situations the 
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Figure 10. 
The form of Iqo as a function of A. 



intermediate region has relatively small size in comparison with the active region (the passive 
region has no size - it can be spread to infinity). Then one can fulfill some further simplifications. 
One can transmit the point w = to the special activity for which 

9 \ w=0 (oo) = - , 

i.e. the half of the centers became the centers of droplets. Then the system of the condensation 
equations after the rescaling of w can be written as 

g = a exp(w) / (z — x) 3 exp(— G)6dx , (69) 

JO 

g(w)dw , (70) 

-oo 

6 = exp(— cli exp(w) / exp(— G)dx) , (71) 
Jo 

where w + is the upper boundary of spectrum and a , ai are some constants. Having rescaled 
z, x one can put a — 1. The condition of the choice for w = gives: 

- In 2 

ai ~ J™exp(-G)dx ' ( ' 

So, there remains only one parameter w+. 
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Let us explicitly extract the intermediate and the active regions. One can see that for two 
different activities W\ and w 2 the following equation is valid 

\n(9(w 1 ,x)) = \n(6(w 2 ,x)) 
exp(wi) exp(-u; 2 ) 

The same is valid for the final values of 9. The value of ln(9(w,x))/ exp(w) is invariant for 
the different activities in one and the same process. So, one can put w = w ++ = 2 4- 3 as the 
boundary between the active region and the intermediate region. In the same manner one can 

separate the region of the passive centers by the boundary w = —w ++ . One can neglect the 

substance in the droplets on the passive centers and get@ 

fW+ 

G ps / g(w)dw . (74) 



One can denote by n QO G + the number of the substance molecules in the droplets formed on 
the active centers. The system of the condensation equations can be rewritten as 



g = a exp(w) / (z — x) 3 exp(— G)dx , 
Jo 

G= g(w)dw + G + , 

J ~oo 

6 = exp(— a\ exp(u>) / exp(— G)dx) . 
Jo 



Now one can write an evident expression for G + 

G + = ^(w + -w ++ )z 3 . (75) 



The system of the condensation equations now is the following one 



g = a exp(w) / (z — x) 3 exp(— G)dx , 
Jo 

/ w ++ Tj^ 
g(w)dw H — w ++ )z 3 

9 = exp(— ai exp(w) / exp(— G)dx) . 
Jo 



The value of w ++ is universal, but the coefficient in the term 77* (w + — w ++ )z 3 /n 00 depends on 
parameters. 

Having defined 

G-= g(w)dw , (76) 

J — 00 

one can propose the following estimate 

G_ < (w ++ -w^)^-z 3 . (77) 



4 The absence of the active and intermediate regions (w+ < w ) means that the condensation occurs in the 

pseudo homogeneous way. 
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Figure 11. 

Behavior of as a function of x for different w. 



When 

w + — w ++ ^> w ++ — w ~ 4 4- 6 

one can approximately get 

G = —(w + - w ++ )z 3 . 

floe 

Then the first equation isn't necessary at all. So, there is no necessity to put ao to 1 and one 
can use the arbitrary scale to cancel another constant. The system looks like 

9 = exp(— ai exp(u>) / exp(— G)dx) , 
Jo 

G = -^-(w + — w ++ )z 3 . 

floe 

Having rescaled x, z, one can put 77* (w+ — = and the system looks like the universal 
expression: 

In 2 f z 

# = exp(-— 5 —- exp(w) / exp(-z 3 )dx) . 

J Q exp{—z 6 )dx Jo 

The behavior of 9 as a function of z for the different values of w is shown in Fig. 11. 

One has to stress that in the general case of quasi continuous spectrum all results concerning 
the intermediate approximate solution can be used. One has only to substitute the sum over 
all sorts of centers by the integral over activities with appropriate normalizing coefficient: 
Ei / dw. 
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Moreover, the procedure of intermediate approximate solution can be directly applied to the 
situation with a mixed quasi continuous spectrum of activities and some discrete 5-like peaks 
for several sorts of heterogeneous centers. 



This remark completes the description of the nucleation for the decay on heterogeneous 
centers with different activities. Unfortunately there are no experimental investigations which 
can be compared with the theory. It is rather natural because earlier there were no formulas 
which could give any physical information and experimental investigations would give only 
some chaotic experimental data. Now one can easily check this theory. This publication can 
be regarded as an invitation to open experimental investigations in this field. One can use 
the same devices as for the ordinary homogeneous nucleation investigations and calculate the 
number of droplets. Experiments on this theory can give information about the properties 
of heterogeneous centers in complex systems. The special interest is to solve some extreme 
problems to produce mainly the droplets on a given sort of centers when there are several sorts 
of centers in the system. With the help of this theory this prolem can be solved. 

References 

[1] Wilson C.T.R. condensation of water vapour in the presence of dust air and other gases - 
Phil.Trans. 1898 V.189A, Nil, p.265 

On the condensation nuclei produced in gases by the action of roentgen rays, uranium 
raeys, ultraviolet rays and other agents - Phil. Trans 1899, V.192A, p. 403 

On comparison efficiency as condensation nuclei of positively and nefatively charged ions. 
Phil. Trans., 1900, V.193A, p.289 

[2] Volmer M., Z.Phys. Chem. vol. 25, p.555 (1929) 

[3] Becker, R. and Doering, W., Ann. Phys. vol. 24, p.749 (1935) 

[4] Zeldovitch, J.B., J.Exp.Theor.Phys. vol.24, p.749 (1942) 

[5] Lothe, J. and Pound, G.M., J.Chem.Phys., vol.36, p.2082 (1962) 

[6] Reiss, H., Katz, J.L., and Cohen, E.R., J.Chem.Phys. vol.48, p.5553 (1968) 

[7] Reiss, H., in A.C.Zettlemeyer (ed.), Nucleation Phenomena, Elsevier New York 1977 

[8] Fisher, M., Physics vol.3 , p.255 (1967) 

[9] D.W.Oxtoby and R.Evans, J.Chem.Phys. vol.89, p.7521 (1988) 

[10] X.C.Zeng, D.W.Oxtoby, J.Chem.Phys, vol 94, 4472 (1991) 

[11] D.W.Oxtoby, R.Talanquer, J.Chem.Phys., vol.100, p.5190 (1994) 



16 



[12] H.Reiss, A.Tabazadech, and J.Talbot, J.Chem. Phys. vol.92, 1266 (1990) 
[13] H.M.Ellerby, C.Weakliem, H.Reiss, J.Chem.Phys. vol.92, p.9209 (1991) 
[14] H.M.Ellerby, H.Reiss, J.Chem.Phys. vol.97, p.5766 (1992) 
[15] C.Weakliem, H.Reiss, J.Chem.Phys. vol.99, p.5374 (1993) 
[16] C.Weakliem, H.Reiss, J.Chem.Phys. vol.101, p.2398 (1994) 
[17] Wakeshima, H., J.Chem.Phys., vol.22, p.1614 (1954) 

[18] Clement C.F., Proceedings of the Royal Society, London A 398, 307 (1985) 

[19] D.R.Warren and J.H.Seinfeld, Journal of Colloid and Interface Science, 105, 136 (1985) 

[20] D.R.Warren and J.H.Seinfeld Aerosol Science and Technology, 3, 135 (1984) 

[21] D.R.Warren and J.H.Seinfeld Aerosol Science and Technology, 4, 31 (1985) 

[22] Barrett J.C., Baldwin T.J., Journal Aerosol Science, 31, 633 (2000) 

[23] H.Reiss, J.Chem.Phys. 19, 482 (1951) 

[24] Barrett J.C., Journal Aerosol Science, 31, 51 (2000) 

[25] V.Kurasov, Phys.Rev. E, vol.49, p.3948 (1994). 

[26] Kuni, F.M., Schekin, A.K., Rusanov, A.I., Colloid.J. USSR vol.44, p.1062 (1982); vol.45, 
p.682 (1983) 

[27] Kurasov V.B., Physica A vol 226, p.117 

[28] Kuni F.M., Novozhilova T.Yu. and Terentiev LA. Letters in Mathematical Physics, vol. 
14, p.161 (1987) 



17 



